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We study the response of the quasi-energy levels in the context of quan- 
tized chaotic systems through the level velocity variance and relate them to 
classical diffusion coefficients using detailed semiclassical analysis. The sys- 
tematic deviations from random matrix theory, assuming independence of 
eigenvectors from eigenvalues, is shown to be connected to classical higher 



£* order time correlations of the chaotic system. We study the standard map 

f— n \ as a specific example, and thus the well known oscillatory behavior of the 

^^ ' diffusion coefficient with respect to the parameter is reflected exactly in the 

o\" 

oscillations of the variance of the level velocities. We study the case of mixed 

phase-space dynamics as well and note a transition in the scaling properties 

O \ of the variance that occurs along with the classical transition to chaos. 
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I. INTRODUCTION 

The quantum spectrum is well known to reflect in several ways classical integrability 
or its lack thereof ]l])@@l- For a completely chaotic, quantized system the energy eigenval- 
ues have characteristic, in fact, universal fluctuation properties that coincide with random 
matrix theory (RMT) universality classes and the eigenfunction components are also dis- 
tributed as Gaussian random variables. However, there are important deviations from this 
dull uniformity imposed by the underlying (asymptotic) deterministic chaos. Classical peri- 
odic orbits, a dense set of measure zero unstable orbits, introduce characteristic deviations 
that are well documented, including the phenomenon of eigenfunction scarring ||. The 
movement of energy levels with the variation of an external parameter, level dynamics, has 
also been studied by several authors with different motivations []5| [3~0|| . It is known that the 



motion of the energy levels as a function of the parameter, now a psuedo-time variable, is 
completely integrable whether the system is itself chaotic or not ||11| . Nevertheless there are 
characteristic features that are introduced by chaos, for instance avoided crossings that may 
be characterized by the second derivative of the energy levels, i.e. the curvatures. 

Here we study level "velocities" , and relate them directly to classical diffusion coefficients; 
although we are using the term velocities, we are not discussing adiabatically changing a 
system, just the slopes of the level curves as a function of a controllable parameter. It has 
been known for some time that these are Gaussian distributed with a variance that has 
been related to a classical "generalized conductance", especially in the context of weakly 
disordered metallic grains. Methods employed were mostly field theoretic and RMT based, 
while numerical simulations of chaotic billiards led to the conjecture that the behavior of 
disordered systems could be extended to chaotic ones as well [|J. 

The variance has a significance beyond setting the scale of the Gaussian distribution of 
velocities. It enters as a normalization required to uncover possible universalities in para- 
metric level correlations. It encodes the system specific characteristics of level motions as 
a function of an external parameter. Level correlations and velocities are experimentally 



accessible, for example in microwave cavities [|12|] or quantum dots, Although, universal 
parametric correlations are not well established experimentally, a recent experiment exploit- 
ing the similarity of elastomechanical wave equations of flexural modes of plates to the 
Schrodinger equation, seems to lend support to it |p~3| . 

In the case when the changing parameters are Aharanov-Bohm flux lines that do not 
lead to any classical dynamical changes, but do lead to important spectral modifications, 
correlation between level velocities were semiclassically considered in 0]. For a treatment 
of Hamiltonian flows see [15|. Recent closely related work, in the context of Hamiltonian 



flows, is also found in ||16|| , where detailed results about the variance of level velocities are 
presented for billiards. 

We make precise the connection between classical diffusion and the variance of the level 
velocities in the simpler context of quantized maps or more generally time periodic systems 
where detailed semiclassical (and classical) analysis is possible. We evaluate the variance for 
the standard map as a function of the external kicking strength and show system specific 
correlations in the form of Bessel function oscillations. Since two-dimensional area preserv- 
ing, or more generally symplectic maps, are Poincare sections of Hamiltonian flows, our 
analysis also reflects upon these systems and is consistent with results derived therein. On 
the other hand, due to the vastly simpler numerical and analytical work involved with maps, 
they lend themselves to more detailed and extensive work. 

We relate our analysis to a semiclassical evaluation of expectation values of generic oper- 
ators in the eigenbasis, as well as touch upon two parameter variations and their correlations. 
The case when the dynamics leads to a mixed phase space is generic and we find a Weyl type 
expansion in h for the variance. The principal contribution in this regime is well predicted 
by a simple classical correlation, which vanishes as the system undergoes a transition to 
chaos. The different scaling behaviors in effective h for mixed and chaotic systems can be 
experimentally observed. We will consider the standard map as an example. Others before 



us have used such systems to study level dynamics P,|TT 



A. The Standard Map and Random Matrix Theory 

Here we define the model studied below and derive the RMT predictions for these. Let 
the classical Hamiltonian have the form 

CO 

H = p 2 /2 - XV (q) £ S(t-n) (1) 

n=— oo 

so that the Floquet operator connecting states just before kicks is given by 

U = exp(-ip 2 /2h) exp(iXV{q)/K). (2) 

The time between kicks is taken to be unity, as there are two independent parameters already 
present, namely A and N. Such systems, known as quantum maps, were first studied in 



|T8| , |T9|| and led to the uncovering of dynamical localization, akin to Anderson localization in 
disordered conductors |21J . We will typically consider the above to be the way the parameter 
of interest (A) enters the problem. 

While this is a map on the plane (for one-degree-of-freedom systems), we consider their 
restriction to the torus [0, l) 2 . This is essential as we have in mind bounded Hamiltonian 
systems and not open scattering ones. Periodic boundary conditions are imposed in both 
p and q directions. We will assume that V(q) is a smooth function on [0, 1) with unit 
periodicity. Denote its average as 

V = f 1 V(q)dq. (3) 

Jo 

Let the quantum map be the N dimensional unitary matrix operator denoted by U. Maps, 

such as the standard map, restricted to a torus are quantized using standard canonical 

quantization [[H|]. Periodic boundary conditions in both canonical variables imposes a finite 

number of states which is the inverse effective Planck constant (h = 1/N). Thus the classical 

limit is approached in the large N limit. Various quantum maps on the torus have been 

studied and form an important part of the literature on quantum chaos due to their inherent 

simplicity P^~p5|. The discrete spectra (N levels) obtained are then analyzed for various 

properties, in particular here the eigenangle velocities are obtained. 



The classical standard map is given by the recursion 

Qi+i = (ft + Pi+i) mod (1) (4) 

Pi+i = (Pi ~ (k/2n) sin(27rgj)) mod (1), 

where i is the discrete time. This is the solution to the Hamiltonian equations of motion for 
the potential V(q) = cos(27rg) and the Hamiltonian in Eq. (|I|). The dynamical variables 
are monitored just before the kicks, and A = k/(27f) 2 . The standard map is of central 
importance as many other maps are locally described by this and the potential may be 
considered to be the first term in the Fourier expansion of more general periodic potentials. 
The parameter k is of principal interest and it controls the degree of chaos in the map, a 
complete transition to ergodicity is attained above values of k « 5, while the last rotational 
KAM torus breaks around k ~ .971 |26|| . 

The quantum map in the discrete position basis is given by |57| 



(n\U\n') = —= exp (irc(n - n') 2 /N) exp j i cos(2vr(n + a)/N) j . (5) 

The parameter to be varied will be the "kicking strength" k, while the phase a will be used 
to avoid exact quantum symmetries, and n, n' = 0, ... ,N — 1. The eigenvalue problem 
of the unitary matrix is written as U\ipj) = exp(—i(f>j)\ipj). The eigenangles <pj are real 
and their variation with the parameter k (level "velocities" ) are given simply by the matrix 
elements 

* = ^s\V\i>i) = ^(^|cos(27rg)|^). (6) 

The 2tt factor is the result of choosing k as the relevant parameter and not k/2n and we 
retain this as this corresponds to the more conventional usage where the last KAM torus 
breaks when the parameter value is just under unity. 

It is then clear that studying level velocities is equivalent to studying expectation values 
of operators in the eigenbasis. Thus if we require (ipi\A\%pi) we would look at the modified 
unitary operator (assuming A is Hermitian) 



U = U exp{-i\A/h) (7) 

where Uq is the quantum system under study. Then the expectation values are simply the 
corresponding level velocities evaluated at A = 0, multiplied by h. If one may identify the 
classical canonical transformation generated by A, we could study a modified classical map, 
as well. However since A = 0, it is the properties of the original classical map that will be 
relevant. The work |28j already discussed the general problem of semiclassical evaluation of 



matrix elements and our following work may be viewed in this context as well. 

From the Gaussian distribution of eigenfunctions for a quantum chaotic system we expect 
the level velocities be similarly distributed. We will concentrate on the variance of these 
velocities, namely the sum: 

We will assume, as is the case with the standard map example, that the average vanishes, 
i.e. V = 0. Later we will generalize to the case of a non- vanishing average, or expectation 
values of operators with non-zero traces. Figure 1 shows a scaled a 2 as a function of the 
parameter k. At about k rs 5 the variance settles down to a near constant, this value 
coincides with the disappearance of major islands of stability in the classical phase space. 
What interests us primarily here however is the clear oscillations that persist as a function 
of k right into regions of large chaos as shown in the inset. 

First we study the value around which the oscillations occur, as this is provided by 
assuming RMT models. Using Eq. @ we get 

nr AT-1 /AT-1 \ 2 



4MT = 4^2 E ( E K^»l V((n + a)/N) 

N 
~ 4^ 



m=0 \n=0 
N N~l /AT-1 

T~ 2 E E K^|n)| 4 (V((n + a)/N)) 2 + (9) 



m=0 \ n=0 



]T \(^ m \n)\ 2 \(i/j m \n')\ 2 V((n + a)/N)V((n' + a)/N) 
The eigenfunctions have been expanded in a basis that diagonalizes the perturbation V, 



which we have taken to be the position basis. Since we assume a zero centered or trace-less 
pertubation 

JV-l 

Y, V{{n + a)/N) = 0. 

n=0 

We use the square of this relation in Eq. (|9J) while replacing eigenfunction components by 
their ensemble averages (denoted by angular brackets) to derive that 

4MT = ^ «l^"»| 4 } ~ (\^m\n)\ 2 \(^ m \n')f)) . (10) 

A crucial step in writing down the above is to assume the independence of the eigenfunction 
components from any specific position eigenvalues. While this is a reasonable statistical 
assumption we will see below that it misses important correlations that are incorporated 
naturally in semiclassical treatments. This is the origin of the non-universality of level 
dynamics, as this implies system dependent correlation effects. The same perturbations (V) 
applied to different chaotic systems will result in different statistical responses, unlike the 
predictions of RMT. 

We use standard results from RMT relevant to the Gaussian Orthogonal Ensemble 
(GOE), which is applicable here as well as the relevant Circular ensembles, [29]. In particular 



(K^H| 4 ) = N{N - 2) ~ jp, (\^m\n)\ 2 \^ m \n')\ 2 ) = N{N + 2) ~ W 
We finally get 

-RMT = ^ (11) 

As a special case for the standard map V 2 = 1/2 and we get Oda/tt = N/Atx 2 . This last 
result explains the value about which the oscillations occur in Fig. 1. This implies that 
the response of the system as measured by the movement of the energy levels is essentially 
the intensity of the perturbation. For chaotic systems then the response is independent 
of the system's detailed dynamical properties. We must also point out that when time 
reversal symmetry is broken the response is half as large. We now turn to the systematic 
oscillations that are not readily predicted by random matrix theory and are manifestly 
system dependent. 



II. SEMICLASSICAL THEORY 



A. The chaotic phase-space 



We first develop in some generality expressions for the variance of the level velocities in 
which semiclassical methods can be easily applied. We write a gaussian smoothed density 



of states BO 



as 



Pm{<S>) = E F M(n) exp(m</>)TrC/ n , (12) 

n=— oo 

where -Fm(^) = exp(— n 2 /2M 2 )/(27r) is introduced to avoid divergences. The exact spiked 
density of states is obtained in the limit M — > oo although almost all levels will be resolved 
at M = N, as the mean level spacing is 2tt/N. The smoothed step function, Nm{4>) is the 
integral of the level density with respect to <fi. We derive then that 



f** ( dN M (<j>) \ 2 _ MNa\k,N) ^ ~ FUn) 
Jo { dk '" 2v^F „f- 



i^ u " 



2 



(13) 



2 V* ' nf^oo n2 

The term n = does not belong in the sum, and it is understood that the first equality is 
an approximation that becomes exact as M — * oo. From this expression it follows that it 
is the long time traces of the propagator, and therefore semiclassically, long periodic orbits 
that are important. 

Another very similar route is through the identity 

JV-l 

Tr(U n V) = Y,(^\V\^)eM-i<t>jn) (14) 



thus implying that 



a\k,N) = ^- 2 (\Tr(U-Vf) n , (15) 



where the angular brackets indicate averaging over time n in the neighborhood of large 
n. We assume as is relevant for chaotic systems that there are no degenerecies. To make 
connections with the standard map above we would take V = cos(27iq). 



Now we make use of the semiclassical approximation of the trace of the propagator as a 
sum over periodic orbits p|j3T[| which is 



Tr (U n ) ~nJ2 Af exp faiNS^ - mv p /2) , (16) 

p 

where A^ 1 ' = l/(2sinh(A p n/2)), and the sum is over periodic orbits of period n which 
labelled by p and have a Lyapunov exponent X p . The actions of these orbits are denoted by 
S^ and are calculated from the generating function of the classical map. The phases v p are 
Maslov like indices and will not concern us here. 

There is also a generalization of the above, which is particularly easy to derive when the 
perturbation is diagonal in the position (or momentum) basis. 

n 

Tr (U n V) ~ J2 A{ p ] ex P faiNS^ - mu p /2) J2 V{tf). (17) 

V j=l 

Here V(x*j) is the value of a phase space representation of the operator V that is evaluated 
along the periodic orbit labelled p and at the point labelled j. An appropriate generalization 



in the energy domain for continuous time systems is found in | 28| . The sum around the 
periodic orbit of the function V is essentially the derivative of the action with the parameter, 
and we may use either the first trace formula in conjunction with Eq. ( |TB| ) or the second 
with Eq. (H). 

We take the second route as we connect with the first subsequently. Taking the modulus 
of the second trace formula gives 

n 

|Tr(trvoi 2 ~£4 n)2 (E^)) 2 + 

V 3=1 

n n 

E 4 n) 4° (E V(#))(E V(qf)) ex P (2mN(S p - S p ,)). (18) 

p=£p' j=l j=l 

As is usual, we have separated the diagonal contribution from the "off-diagonal", which 
corresponds to distinct pairs of orbits, with distinct actions. We have also assumed for 
simplicity, as is the case with the specific parameter variation chosen above in the standard 
map, that V(x) is only position dependent; this does not alter the results below. We have 
also included the phases into the actions. 

9 



Since we expect that long periodic orbits are important, the diagonal approximation, 
which relies on random phases may be violated due to subtle correlations among their 
actions. The time at which we may expect action differences of the order of H is the so 
called log-time, or Ehrenfest time. We argue that action differences are of the order of the 
orbit separation, and since areas of the order of % (for two-dimensional maps) would be 
populated with multiple periodic orbits beyond the log-time, their action differences would 
also be comparable with %. However long periodic orbit actions are randomly distributed 
and will acquire correlations only around the Heisenberg time. At this time the off-diagonal 
terms will dominate the sum, as happens if we simply consider < |Tr([/ n )| 2 >„, which is 
asymptotically N, while the diagonal term is linearly increasing in time. 

However the off-diagonal term vanishes due to the sums of V(q) over very long periodic 
orbits. We may write for two distinct orbits p and p' after assuming uniform measure and 
replacing time averages over the periodic orbit, by the phase-space average that 

n n 

(En^))(En<))~^ 2 (19) 

3=1 3=1 

From our initial assumption that V = 0, the off-diagonal term vanishes. The diagonal term 
is non-vanishing as we once again treat periodic orbits as ordinary long chaotic trajectories 
and derive that 

1 n n 

-(£ V(qf)) 2 = D(k) = C(0) + 2 £ C{1) (20) 

U 3=1 1=1 

where the time correlations are replaced by classical phase space averages due to ergodicity. 

C(l) = (V(q )V( qi )) = ^J A dq dp V(q )V(f l (q ,p )), (21) 

where f l (qo,Po) — qi, f l is the integrated dynamics in time /, and A denotes both the 
phase space and its area (in the cases considered this is unity). We assume that these exist, 
and are decreasing with /, typically exponentially for chaotic systems and that a few terms 
may be sufficient. This is not established in generality and complications may arise due 
to marginally stable orbits leading to non-exponential behaviors. For the standard map, 

10 



coefficients upto C(2) are dominant and sufficient to see the essential behaviour. We have 
dropped the index p as now we will treat such long periodic orbits as generic non-periodic 
orbits. Indeed by using the ergodic theorem we have already abandoned any particularities 
that may arise due to the orbit being periodic. Later, we remark on a case when we may 
not neglect off-diagonal terms. 

The alternative route is to take the derivative of the first trace formula and use Eq. (|13|). 
Again we neglect the off-diagonal terms for reasons given above. We will assume that the 
derivatives of the actions with the parameter ("action velocities") , for a given period or 
period interval, are such that their average is zero while their variance is proportional to 
the time period. This assumption is equivalent to the vanishing of the phase space average 
of V(q) and the presence of ergodicity. This was noted for general Hamiltonian systems in 
|32|1 , and we will see below in the context of maps how this simply arises. We replace then 



for each time n the individual action velocities (squared) by the variance, 

(Pf)V DW "' (22) 

where the angular brackets indicate the average over periodic orbits of period n. 

In either approach, the uniformity principle |33[ is applied in the form that there are e hn /n 



orbits each with a Lyapunov exponent approximately A per unit time. Then |A^n)| 2 ~ e Xn , 
and assuming near equality of the topological entropy h and A, we derive from Eq. (|18D that 



|Tr(fW)| 2 ~ gD(k). (23) 

Similarly from the other approach 



n 



4rTrU n 
dk 



' N 2 



9^- 2 D(k). (24) 



The tilde sign in the above equations implies that the L.H.S can be expected to be the 
R.H.S in an average sense. The spread in time n will also reflect the spread in the average 
action velocity diffusion coefficent D(k) with period. Results not shown here indicate that in 
the chaotic regime this is an exponential distribution. The factor g inserted above is due the 

11 



fact that symmetries can impose distinct orbits to have identical actions. This factor must 
be determined from classical and quantal symmetries, and includes phase-space symmetries. 
We finally get then from either approach the response in terms of the variance of the level 
velocities: 

a 2 (k,N) = g JL D {k). (25) 

Thus the variance of the level velocities is proportional to a classical diffusion coefficient 

that determines the diffusion of action velocities of periodic orbits. More explicit expressions 

for this coefficient are now derived. Area preserving maps, such as the standard map, have 

a generating function L(qi+i,qi] k) from which the map may be derived as dL/dqi = —pi 

and dL/dqi+i = p-i + \ (I. C Percival in [^J). The total action of a periodic orbit is equal to 

S^ = ^2iL(qi + i,qi; k), where the sum is over the n periodic points of the orbit p. Thus 

we derive, after assuming that the orbit is not at a point of bifurcation, that for a periodic 

orbit: 

dS^ dS^ 
dk ~ dk ■ [Zb) 

We have not used the partial derivative sign in defining the level velocities although we 

assume that only one parameter is varied. This is due to the subsequent fact that when the 

classical action derivative is written, it is a total derivative, in as much as the periodic orbit 

itself changes with the parameter. These two, however, are shown to be equal in the case of 

periodic orbits. 

The variance is given by 

((^Ir) ) = (.E i cos ( 2 ^)cos(27rg i )^ ~ |(1 + 2J 2 (A;)) (27) 

In the equality the sum is over different times along a given periodic orbit and then averaged 
over all periodic orbits of period n, while the approximation arises from a replacement of the 
average by the usual ensemble average and retaining up to the second order time correlation. 
More precisely C(0) = 1/2, C(l) = 0, C(2) = J 2 (k)/2, where J 2 (k) is a Bessel function. 
These are derived from: 
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I n pi pi 

C(0) = -Y cos 2 (2vrg,) ~ / / cos 2 (2vrg) dqdp = 1/2, 
n ~^ Jo Jo 



1 " 
C(l) = - Y, cos(27rgj) cos(27rg i+ i) 



n r=i 






/" f cos(2nq)cos(2ir (q + p-(k/2ir) sm(2irq)))dq dp = 0, (28) 

1 n 
C(2) = - V cos(27rgj_i) cos(27rg i+ i) 

~ / / cos(27r(g-p))cos(27r(g + p- (A;/2tt) sin(27rg))) rfgrfp = J 2 (k)/2. 
Jo Jo 

The symmetry factor is g = 2, for the standard map, if we assume time-reversal invariance 
alone. This is the case for the data presented in Fig. 1 as we have intentionally broken the 
phase space symmetry in the quantum system by assuming a = .35, (generic value) rather 
than a = .5, which will lead to twice the variance. Periodic orbits are either self-symmetric 
or more generically have symmetric partners with identical actions. 

The Bessel functions are characteristic of the diffusion coefficient in the standard map, 
the above simple derivation was proposed in the context of deterministic diffusion in p| 



The linear time dependence is then a consequence of ergodicity, whereby time averages are 
replaced by phase space averages, and the coefficient is also easy to find and is the scaling of 
the parameter introduced to uncover possible universalities in level dynamics. These then 
are the oscillations observed in the figure. The bold line in the inset is 2D (k) = (l + 2J2(/c)). 
The significant deviation around the first minimum in the inset (k ~ 6.5) from theory could 
be due to the presence of small stable islands, which are the accelerator modes and are 
known to lead to anomalous transport in the standard map (B. V. Chrikov in ||). 

B. The mixed phase-space 

The regime where there is a mixed phase-space consisting of large stable regions is generic, 
and in this case the analysis above fails: the assumptions about the trace formula and the 
uniformity principle operate only under conditions of complete hyperbolicity. While in the 
completely chaotic regime the variance scales as N, in the mixed phase space regime it 
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(principally) scales as iV 2 . This relates to the large hump in Fig. 1, to which we now turn 
our attention. One way of relating the variance to classical quantities is to recognize that 

a 2 (k,N) = — (Tr(U- n VU n V)) , (29) 

47T \ In 

where the average is once more the time average. Thus the variance of the quantum level 
velocities is directly the time average of operator auto-correlations. We consider the case 
when the average level velocity is zero, as the generalization is evident. Replacing opera- 
tors by the corresponding classical observables, we expect to get the variance in the mixed 
phase regime. This is particularly successful as we are dealing with averages over the entire 
quantum spectrum. Thus we replace the trace operation divided by iV with the classical 
phase-space average to get 

a 2 (k,N) cl = —j o dp j o dq V(q ) (V(? (n) («b,Pd))) B , (30) 

where q^ n \qo,Po) is the position after n iterations starting from the initial condition (qo,Po)- 
For the case in Fig. 1, V(q) = cos(2irq) and Fig. 2 compares in the mixed phase space 
regime the exact quantum calculation with a purely classical simulation corresponding to 
a 2 (k,N) c i. We see that a simple classical simulation reproduces the curve extremely well, 
including the secondary hump, till around k rs 2tt. It is quite remarkable that the classical 
curve continues to pick out the initial bessel function oscillations in the deeply chaotic regime. 
In the figure the time average is done over an ensemble for the first hundred iterations, and 
the oscillations indicate short time correlations that will strictly disappear with increasing 
time. 

The transition to classical chaos is accompanied then by a transition of the variance of 
the level velocities from a quadratic to a linear iV dependence. Based on this observation 
we may write a general expression for the variance of the level velocities as a Weyl series 
with principal terms 

a 2 (k,N) = Cl (k)N + c 2 (k)N 2 , (31) 
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where c\(k) and 02(h) are system dependent and we have given above their expressions 
assuming only one of them is appreciable. Note that we have not evaluated Ci(k) in the 
mixed phase-space regime, and that this will not in general vanish. On the other hand we 
expect that c^ik) vanishes as a classical transition to chaos occurs. This is illustrated in 
Fig. 3, where c 2 (fc) is evaluated based on a best fitting curve using five N values, equally 
spaced, between 100 and 500. The curve is fitted by assuming a third-order polynomial in 
N for which the co-efficient of N 3 returned by the fit was always of the order of 10~ 6 or less. 
In general the RMT result derived earlier Eq. (£Ll]) will be correct under the assump- 
tion that C(l) = C(0)So t i, implying delta correlated processes. Thus the departures from 
universality is related to higher order time correlations. The response of the system is not 
only dependent on the strength of the perturbation, but also on the dynamical correlations 
inherent to the system. 

C. Generalizations 

We remark now on the general case V 7^ 0. This implies an overall drift to the energy 
levels due to changing phase space volumes. Using Eq. fll9|) and adding and subtracting 
nV from the diagonal part of Eq. flllf ), we get after using < |Tr (U n )\ 2 > n = N, 

a 2 (k,N) = ^ g (D(k) - V 2 ) . (32) 

The large chaos limit of this is 

a\oo,N) = ^- 2 g(W-V 2 ), (33) 

and is the RMT result. 

Variations of two independent parameters is an important problem, considering that 
many novel effects, including geometric phases may be observed. Here we will consider, in 
a generalization of the above, correlations between independent parameter variations. We 
assume that the Level velocities are given by the matrix elements: 
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i = 1,2, and Aj are two independent parameters while VJ are two Hermitian operators. Thus 
the correlations considered below are also correlations between diagonal matrix elements of 
two arbitrary Hermitian operators. 
We derive then that 

Using methods as outlined above the correlation function is semiclassically evaluated and 
we get 



^ Ai ' A2 ) = r^((i?if ) "(if) (if) I - ( D(Al,A2) VlV2 \ (36 ) 

(<ri<7 2 ) \\dX 1 dX 2 / j \dX 1 / j \d\ 2 / j J .^ Dl _v[ 2 )(D 2 -V2 2 ) 

The function D is a generalization of Eq. (^) involving the dynamical correlation between 
the functions V\ and V 2 . 

D(Ai,A 2 ) = -j ( / ^i(g ,Po)^2(go,Po)rfgorfpo + 2 J2 J v i{<lo,Po)V2{qi,Pi)dq dpo\ , (37) 

the dynamical variables after a time / integrated from (qo,po) is denoted (qi,pi). While Z?i,2 
refer to the correlations appropriate to them individually and defined earlier in Eq. (p0|) . 
We remark that our derivations have assumed time-reversal symmetry, and that the factor g 
is responsible also for phase-space symmetries. It may be generalized to include the factors 
that come due to breakdown of time-reversal, or inclusion of spin. 

We finally consider the situation where not only the averages V i vanish, but also there 
is no tangible correlation between them, i.e., DiX^X-i) = 0. Then the semiclassical ex- 
pressions above give zero and are incapable of capturing the small albeit finite and rapid 
oscillations (with parameter). This limitation is of course evident all along, including Fig. 
1, where the Bessel function captures the low frequency oscillations only. 

In Fig. 4, is one such example, where we have considered V% = cos(2irq) and V 2 = 
cos(47rg). The parameter X 2 is set as zero so that the relevant classical system is still the 
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standard map, while the other parameter is k above. The correlation is seen as a function of 
this last parameter, the average actually vanishing. We see that this measure too captures 
the transition from mixed phase space to chaotic phase space, but that after the transition 
there are only extremely rapid oscillations about zero, although there are quite frequently 
fairly large correlations. In fact the frequency of the oscillations are so rapid that they seem 
to have self-similar properties as a random fractal. We may estimate the order of magnitude 
of the frequency if we assume that these arise from the off-diagonal part of the semiclassical 
sums. The magnitude of the parameter change needed to change a typical orbit action by h 
or (1/N) is needed. From the fact that action changes have a variance proportional to the 
period, we get \AS\ ~ \/n|AA|, and therefore |AA| ~ iV~ 3 / 2 , if we take as the period n = N, 
which is the Heisenberg time and represents the time by which the spectrum is practically 
resolved. 

In conclusion then, we have studied variances of level velocities and their generalizations 
in the chaotic as well as mixed phase space regimes. Noting that the transition to chaos 
is perfectly reflected in this measure, we derived detailed formulae for them, in terms of 
classical correlation coefficients and illustrated this with the help of the standard map. The 
mixed phase space regime was surprisingly well captured by a simple classical estimate. 
The observations of oscillations or variations in the level velocity variances due to classical 
correlations, as well as using them to distinguish mixed from chaotic phase space are both 
experimentally accessible. 

The possibility of using the level velocity in conjunction with wavefunction intensities in a 
measure of phase-space localization has been proposed f35fl , and it is hoped that this detailed 



understanding of the level velocities will help in this as well. In particular this measure is 
a special case of the correlation between two operators discussed above, with an important 
complication being that the Wigner transform of the relevant operator, a projector in phase 
space, varies over scales of order h. We have also noted that the RMT results, after assuming 
independence of eigenvalues and eigenfunctions is capable of predicting the level velocities 
only in the limit of extremely large chaos, or equivalently ignoring all higher order time 
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correlations other than the zeroth. 
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Figure Captions 

Figure 1 Scaled variance as a function of the parameter k, N=300, a=0.35. The inset 
shows a part of the plot magnified, the points are numerical data while the smooth line is 
the twice the diffusion coefficient. 

Figure 2 Same as Fig. 1, comparing, mainly in the mixed phase space regime, the 
exact variance (dots) with the classical estimate (the line). The classical estimate is after 
averaging over a hundred time steps. 

Figure 3 The coefficient 02(h) as a function of the parameter; note that it reflects the 
classical transition to chaos. 

Figure 4 The quantum correlation between the matrix elements of the two operators 
cos(27rg) and cos(47rg) whose classical correlation vanishes in the chaotic regime. 
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